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The absolute dose difference of the iso-center position at different times was used to evaluate the stability 
of medical linear accelerator (LINAC). In order to eliminate the drawback of point dose without statistical 
significance and to investigate the output stability of segment with small weight related to intensity-modulated 
radiotherapy (IMRT), a modified method of evaluating the stability of LINAC was proposed based on planar dose 
combined with gamma method. With a commercial available two-dimensional ionization chamber, a set of planar 
doses with dose gradient from 2 cGy to 100 cGy were obtained. Then the gamma method was adopted to analyze 
the dose difference between the reference data and the evaluated data at each dose level. The results showed that 
the improved method based on planar dose for analyzing the output stability of the medical LINAC was feasible 
and efficient, and suggested that the reverse optimization should be aborted in clinical when the segment weight 


related to IMRT was under 10 MU. 
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I. INTRODUCTION 


Output stability monitoring of medical LINAC is an im- 
portant aspect of radiotherapy quality assurance. The dai- 
ly/weekly check for output stability of medical LINAC cur- 
rently in clinical practice was carried out according to the pro- 
cedures of the International Atomic Energy Agency (IAEA) 
report No. 227 [1]. The dose difference of the iso-center point 
at different time was adopted to evaluate the output stability. 
However, with the clinical implementation of new dose mea- 
surement equipment, point dose showed the drawback that it 
could not reflect the planar dose output stability of the LINAC, 
so it had no statistical significance. This was definitely the 
case of Intensity-modulated Radiotherapy (IMRT), in which 
the planar dose evaluation for segment with small weight is 
of great importance. The gamma method introduced by Low 
Daniel A et al [2] was the most accepted approach and had 
been widely used by many groups [3—5]. In this paper, a mod- 
ified evaluation method of the stability of medical LINAC was 
proposed. The gamma method was adopted to do planar dose 
verification and the output stability of LINAC was evaluated. 
This research was a part of Advanced/Accurate Radiotherapy 
System (ARTS) [6-12], in precision radiation treatment plan- 
ning and quality assurance system project, developed by FDS 
Team [13-19]. 
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II. MATERIALS AND METHODS 
A. LINAC and MatriXX 


The measurements were conducted using photon beam from 
XHA600D (SHINVA, China)(see Fig. 1(a)). The LINAC could 
produce 6 MV X-rays. Using the three-dimensional (3D) Blue 
Water Phantom (IBA, Germany), the maximum dose depth was 
detected to be 1.5 cm under water. The I’mRT MatriXX (IBA, 
Germany)(Fig. 1(b)) device consists of a two-dimensional (2D) 
array of ionization chambers. There are 1020 vented parallel 
plate ionization chambers on the array detector, arranged in 32 
x 32 grid. The center-to-center distance of chamber is 7.62 
mm, with an active area of 24 cm x 24 cm. The MatriXX was 
positioned using the LINAC field light. The beam buildup of 
MatriXX was about 3 mm. On the MatriXX surface, 30 cm x 
30 cm RW3 Solid Water (Sun Nuclear, America) layers served 
as beam buildup with a thickness of 1.2 cm. The LINAC 
gantry angle was set to be 0 degree, and the source-to-surface 
distance (S.S D) was 100 cm. 


B. Calibration and measurement 
1. Calibration 


As recommended by international standards[1], under the 
radiation of SSD = 100 cm and the field size is 10 cm x 10 
cm, the temperature T = 24.6°C, the pressure P = 999.7Pa, 
the LINAC dosimeter should be modulated the output of 100 
MU (MU: Monitoring Units) proportional to 100 cGy using 
the cylindrical ionization chamber with 0.125 cm? air volume 
(PTW TM 31010 S/N03422) at the maximum dose depth 
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(1.5 cm under the solid water for 6 MV photon). 


The MatriXX was put under the measurement condition 
for 30 min in order to do the background correction. The 
solid water with 1.2 cm thickness was added to the surface 
(the device had a beam buildup of 3 mm). At each MU level, 
the average value was calculated for the center four pixels. 
The cylindrical ionization chamber was used to measure the 
absolute dose. After that the calibration scaling factor k was 
calculated and saved in the software of Omni I'm RT. 


2. Data measurement 


Under the radiation measurement condition (SSD = 100 
cm, field size: 10 cm x 10 cm), after LINAC warming up for 
30 minutes (marked as period 1), different weight were irradi- 
ated on 1.2 cm solid water with 1 MU, 2 MU, 3 MU, 5 MU, 
10 MU, 20 MU, 40 MU, 50 MU, 100 MU, respectively. The 
planar absolute/relative doses were obtained with MatriXX. 
For each MU level more than 10 planar doses were collected 
and the average planar doses were calculated, which was used 
as the reference data for gamma method. After keeping power 
on for 8 hours (marked as period 2), the same protocols were 
adopted and the average planar dose was used as the evaluated 
data. 


With SSD = 100 cm, 28 cm x 28 cm field size and 100 
MU was preset, the active area of 24 cm x 24 cm (maxi- 
mum field) was covered to investigate the beam uniformity 
of LINAC. The flatness and symmetry were key performance 
indicators of beam dose distribution. According to the Inter- 
national Electrotechnical Commission (IEC) standards, the 
flatness and symmetry of beam field should be within +3%. 
10 different X axis position and 10 different Y axis position 
were chosen to verify the flatness and symmetry of the maxi- 
mum field. Table 1 showed that the flatness and symmetry of 
all positions were within +3% and dose distribution was rel- 
atively uniform and there was no obvious noise point, which 
indicated that measured data obtained from MatriXX were 
reliable and could be used for dose verification. 


Table 1. The flatness and symmetry of different X axis and Y axis 
of the maximum field. 


FLATNESS SYMMETRY 

X axis (%) Y axis(%) | X axis(%) Y axis (96) 
1 102.70 101.97 100.78 100.51 
2 102.72 101.99 100.78 100.49 
3 102.72 101.99 100.84 100.54 
5 102.72 102.00 100.81 100.51 
6 102.72 102.00 100.80 100.52 
7 102.69 101.97 100.77 100.46 
8 102.71 101.95 100.80 100.49 
9 102.69 101.97 100.79 100.48 
10 102.70 101.92 100.77 100.51 
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Fig. 1. (Color online) LINAC(a) and MatriXX (b) 
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Fig. 2. Gamma method realization flowchart with VC++ 


C. Gamma (y) method 


The gamma method was mainly used to do the dose verifi- 
cation and analyze the difference between the reference data 
and the evaluated data. Expressed in formula, it could be read 


Ar(r,r').2 AD(r,r').2 
x ew [ert AMM o 


oro 


AD : Dose difference between the reference data and evalu- 
ated data. 

Ar : Spatial distance of the reference data and evaluated data. 
do: Distance to agreement tolerance, usually preset 3 mm/2 
mm in clinical practice; 

9 Do: Dose difference tolerance, usually preset 596/396 in clin- 
ical practice. 

The pass-fail criteria for gamma index for certain point was 
as following: 

when »(r) < 1, pass, 

when »(r) > 1, fail. 

Then y(r) values of all evaluated data consisted of thedis- 
tribution. If the value was less than or equal to 1, the point 
was considered to pass the verification. The pass rate was 
obtained through the pass points number divided by the total 
number of the evaluated data. If the pass rate was more than 
95%, it represented the good agreement of two planar dose 
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Table 2. The readout of cylindrical ionization chamber at different time 
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MU’ Data! ”(cGy) Data2 *(cGy) Error 1*(96) Error 2 *(96) Aerror (96) 
MAX MIN MEAN MAX MIN MEAN 
1 1.536 1.406 1.459 1.458 1.388 1.418 45.8 41.8 5.10 
2 2.582 2.544 2.562 2.548 2.414 2.449 28.1 22.4 5.65 
5 5.572 5.526 5.542 5.501 5.397 5.415 10.84 8.30 2.54 
10 10.64 10.68 10.658 10.65 10.64 10.642 6.58 6.42 0.16 
20 20.74 20.68 20.702 20.78 20.73 20.756 3.51 3.78 0.27 
30 31.03 30.78 30.80 31.04 30.76 30.894 2.67 2.98 0.30 
50 50.92 50.84 50.89 51 50.93 50.97 1.78 1.94 0.16 
80 81.16 80.04 80.9 81.06 81.19 81.1 1.125 1.375 0.25 
100 101.3 101 101.16 101 101 101 1.16 1.0 0.16 
à monitor units, output weight for medical LINAC; 
P 10 times readout of “pinpoint” chamber (period 1); 
* 10 times readout of “pinpoint” chamber (period 2); 
d relative error between the mean measure and MU (period 1); 
* relative error between the mean measure and MU (period 2); 
f Aerror, relative error of mean dose between two periods. Aerror = 10096*(Dataly gAN -Data2y gAN)/MU. 
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Fig. 3. (Color online) Output dose fluctuation of (a) 1 MU, (b) 2 MU, (c) 5 MU at period 1 
t : calculate the value f: h point in the evaluated 
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data. 


MU Pass rate (%) Pass/Fail 
1 20.83 Fail 
2 39.06 Fail 
5 64.37 Fail 
10 93.03 Pass 
20 95.84 Pass 
30 97.59 Pass 
50 95.32 Pass 
80 95.97 Pass 
100 100 Pass 


Step 4: record the pixel number with the value less than or 


equal to 1. 


Step 5: calculate the pass rate of the evaluated data. 


III. 


RESULTS 


A. LINAC output stability for different MU level 


distributions [20, 21]. Since the resolution of the evaluated 
data was limited and gamma method required no less than 1 
mm data resolution, we used the tool of Microsoft VC++ to 
realize the dose interpolation and dose comparison according 
to the formula (1) in this paper. The tolerance of 396/35 mm 
was chosen. The flowchart of the program realization was just 
as Fig. 2. 


Step 1: input the data exported from Omni I’m RT with the 
format *.txt". 


Step 2: interpolate the evaluated data making the data reso- 
lution from 32x32 to 320x 320. 


Under the radiation condition of SSD = 100 cm, 10 
cmx 10 cm field size, temperature T = 24.6°C, pressure 
P = 999.7 Pa, calibration factor ND = 1.026, the 100 MU 
output of LINAC was modulated proportional to 100 cGy. 
The cylindrical ionization chamber (*pinpoint" chamber) was 
adopted to measure the output dose at the position under 1.5 
cm of solid water phantom. Table 2 showed the mean readout 
of 10 times and the output difference for each MU level. 

After MatriXX calibrated, the planar relative dose signal 
could be collected using the movie mode (collected time = 
500 ms) and exported as document with the format ".txt". The 
scaling factor k gained from absolute dose calibration was 
multiplied by the relative dose distribution, then we could get 
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Table 4. The relative dose distribution in 10 cmx 10 cm open field 
X -4.191 -3.429 -2.667 -1.905 -1.143 -0.381 0.38 1.142 1.904 2.666 3.428 4.19 
Y 
-4.191 23 24 23 23 23 23 23 23 23 23 14 13 
-3.429 24 24 24 23 23 24 23 23 24 24 23 23 
-2.667 24 24 24 23 23 24 23 23 24 24 23 23 
-1.905 24 24 24 23 24 24 23 23 24 24 24 23 
-1.143 24 24 24 23 24 24 23 23 24 24 23 23 
-0.381 24 24 23 23 23 24 23 23 24 24 24 23 
0.38 23 23 24 24 23 23 24 23 23 24 24 23 
1.142 23 23 24 24 23 23 24 24 23 24 23 22 
1.904 23 23 24 24 23 23 23 23 23 24 23 22 
2.666 23 23 24 24 23 23 24 23 23 24 23 22 
3.428 23 23 24 24 23 23 23 23 23 24 23 22 
4.19 22 23 23 23 23 23 23 23 23 23 23 22 


the planar absolute dose. The planar dose obtained at period 1 
was input into the gamma method program as the reference 
data while the data at period 2 as the Evaluated Data. Later 
the pass rate was calculated. Table 3 showed the absolute 
dose comparison results of two different periods for each MU 
gradient, and pass means that the dose distribution was similar 
or thought to be the same, while fail representing there was 
obvious difference between two dose distributions. The 3%/3 
mm tolerance was chosen in this study. 


B. Segment with small MU 


In clinical practice, segments were adopted for IMRT and 
every fraction dose was usually optimized to be several seg- 
ments with different weight. In this paper the segment with 
small MU output stability of LINAC was investigated. As in 
the Table 4, although the measured data of relative dose for 2 
MU showed good uniformity, when multiplied by the scaling 
factor k (0.29), the absolute dose distribution was obtained, 
then the relative error of output dose was calculated. 

[((24 x 0.29) — 2]/2 x 100% = 248% 

It was obvious that segment with small MU was fluctuated 
from the calibrated dose greatly. However, this was just the 
situation of small MU segment. From the results in Table 
3 it could be concluded that when the segment with small 
weight the dose obtained failed the dose comparison, for ex- 
ample, 1 MU, 2 MU and 5 MU failed the comparison. When 
the segment weight was above 10 MU, the dose error was 
accepted. 

The readout data of cylindrical ionization chamber was used 
to verify the result, as in the Fig. 3, X axis represented the 
measure times and Y axis represented readout data. Period 1 
was the time after LINAC warming up. It is obviously that the 
output weight fluctuated severely under 10 MU, and deviated 
from the calibrated data. The same trend of instability was 
also found at period 2. It inferred that small MU related to 
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IV. DISCUSSION AND CONCLUSION 
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